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ABSTRACT 

Context. Molecular clouds near the Hii regions tend to harbor more luminous protostars. 

Aims. Our aim in this paper is to investigate whether or not radiation-driven implosion mechanism enhances luminosity of protostars near 
regions of high-ionizing fluxes. 

Methods. We performed numerical simulations to model collapse of cores exposed to UV radiation from O stars. We investigated dependence 
of mass loss rates on the initial density profiles of cores and variation of UV fluxes. We derived simple analytic estimates of accretion rates and 
final masses of protostars. 

Results. Radiation-driven implosion mechanism can increase accretion rates of protostars by 1-2 orders of magnitude. On the other hand, mass 
loss due to photo-evaporation is not large enough to have a significant impact on the luminosity. The increase of accretion rate makes luminosity 
1-2 orders higher than those of protostars that form without external triggering. 

Conclusions. Radiation-driven implosion can help explain the observed higher luminosity of protostars in molecular clouds near H ii regions. 
Key words. Stars: formation - ISM: HII regions - Method: numerical 



1. Introduction 

Massive stars play a very important role in star formation his- 
tory in giant molecular clouds. Strong radiation from massive 
stars impacts on surrounding environment, and may trigger star 
formation nearby. Two scenarios of triggered star formation 
near Hii regions h ave been suggested. One of them is "col- 
lect a nd collapse" dElmegreen & Ladalll977t IWhitworth et al 



19941) . As the H ii region expands, it pushes ambient gas into a 
shell. Subsequent fragmentation and collapse of the swept-up 
shell may form the next generation of stars. For example, a core 
within one fragment of the molecular ring surr ounding Galactic 



Hi: region Shi 04 contains a young cluster dPeharveng et al 



12003) . Ra diation-driven implosion (hereafter RDI) is the other 
scenario dBertoldil [l989l iBertoldi & McKeelll990h . When an 
expanding Hii region engulfs a cloud core, the strong UV 
radiation will compress the core. Bright-rimmed clouds of 
cometary shapes found at the edges of relatively old Hii re- 
gions, are explained bv iLefloch & Lazareffi (119941) with a RDI 
model. To understand star formation in giant molecular clouds, 
we study effects of the RDI. 

Observational results indicate t hat the Hii regions can 



strongly affect star formation nearbv. lDobashi et al.l d2001h in- 
vestigated the maximum luminosity of protostars as a function 



of the parent cloud mass. They found that molecular clouds 
near H ii regions contain more luminous protostars than others. 
Moreover, Sugitani et al. ( 1989) showed that the ratios of lu- 
minosity of protostar to core mass of bright-rimmed globules 
are much higher than those in dark globules. These results in- 
dicate that the RDI affects the luminosity of protostars near H n 
region. 

Recent studies suggest that triggered star formation by ex- 
ternal compres sion can increase accretion ra tes and luminosity 
of protostars. iMotovama & Yoshidal d2003b investigated col- 
lapses of cores compressed by external shock waves and found 
that accret ion rates were enhan ced. Some observations support 
this result. iBelloche et"al] d200 6) observed the Class protostar 
IRAS 4A in star formation cloud NGC1333, whose star forma- 
tion is suspected to be tr iggered by powerful molecular out- 
flows ( Sandell & Kneel2001i) . They showed a numerical model 
of collapse triggered by external pressure, which reproduces 
the observed density and velocity profiles. Moreover, they de- 
duced a high accretion rate of (0.7-2)x lO"'* MQyr"'. Their re- 
sult indicates external compression increases the accretion rate 
of a protostar, which produces high accretion luminosity. In this 
paper, we attribute the origin of higher luminosity in protostars 
to star formation activities induced by the RDI. 
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Some studies h a ve bee n made on RDI. iBertoldi (Il989h and 
Bertoldi & McKed (Il990l) developed an approximate analyti- 
cal solution f or the evolution of a cloud exposed to ionizing 
UV radiation. ILefloch & LazareffI ( Il994 investigated the RDI 
with numerical simulations. However, their calculations did not 
include self-gravity of the gas. Rec ently, S PH simulations in- 
cluding self-gravity have been done. iKessel -Deynet & Burkert 
(l2Q03h mvestigated the effects of initial density per turbations 
on the dynamics of ionizing clouds, and^ Miao et al. (1200 6) fo- 
cused on cloud morphology. Despite of many studies on RDI, 
little attention has been given to the effects of RDI on accretion 
rates. 

In this paper, we explore if the RDI can enhance accretion 
rates by orders of magnitude compared with those of proto- 
stars that form without external trigger. RDI is expected to en- 
hance the accretion rates as a result of strong compression, and 
it would increase luminosity of protostars. On the other hand, 
photo-evaporation of the parent core due to RDI decreases the 
mass of protostar, subsequently decrease the luminosity of the 
new protostar We investigate which effect dominates. 

The paper is organized as follows. In Section 2, we describe 
our numerical method. In Section 3, we present our results. In 
Section 4, we discuss luminosities of protostars based on our 
numerical results. In Section 5, we summarize our main con- 
clusions. 



2. Numerical Method 

We simulate evolution of cores exposed to diffuse UV radia- 
tion. We assumed that a core of neutral gas is immersed within 
an H 11 region and is exposed to diffuse radiation from the sur- 
rounding ionized gas. In an actual Hii region, the core is not 
only exposed to diffuse radiation but also to direct radiation 
from the ionizing star. The diffuse flux of Lyman continuum 
photo ns is ~ 15% of the direct flux of L yman continuum pho- 
tons JCanto et al.l Il998t IWhitworth & Zinnecked 120041) . The 
spherical core has an effective cross-section to diffuse radiation 
4nr^ and to direct UV radiation nr^, respectively. The ratio of 
total number of diffuse photons to total number of direct pho- 
tons is ~ 0.6. Our assumption underestimates the total amount 
of ionizing UV photons by a factor of ~ 2.6. Since the purpose 
of this paper is to investigate the effects of RDI on the accre- 
tion rate of a protostar, we neglected direct radiation and con- 
sidered only diffuse radiation for simplicity. Effects of direct 
radiation and morphology of the core, are beyond the purpose 
of the present paper, and will be discussed in our next paper 
using two dimensional simulations. 

The ionizing UV photons enter from outer boundary of the 
computational domain. The number flux of these UV photons 
is expressed as 



Fi = 0.15 X 



(1) 



where A^^ is the total flux of ionizing UV photons from the star, 
and d is distance from the ionizing star to the core. We assume 
that the ionizing star is an 08V or an 04V star, whose lo R Nl 
is 48.87 s ' and 49.70 s ^ respectively dVacca etalJll996h . 



We assumed that the sound speed is a function of the ion- 
ization fraction x. The hot ionized gas has a temperature of 
~ 10"* K in an Hn region. To imitate condition of the hot 
ionized gas (x = 1), sound speed of the ionized gas c, was 
set to 13kms"'. The neutral gas (x = 0) has two states: 
warm gas, if the density is lower than a critical density pcrit - 
1.67 X lO"^-' (10 ihh) g cm"-', where nifj is the mass of hydro- 
gen atom, and cold gas, if the density is higher than the critical 
density pcrii- We introduce this artificial condition for the ini- 
tial pressure balance between the core and the rarefied warm 
gas surrounding it. The rarefied warm gas is ionized in a very 
short time and will hardly affect the time evolution of the core. 
The sound speed of neutral gas c„ is expressed as 



Ccokl - 



- 10 kms ' if p < p, 
0.2 km s"' if p > p, 



cnt 
'crit 



(2) 



where c„,n„„ and Ccou are the sound speeds of warm gas and 
cold gas, respectively. Sound speed of the partially ionized gas 
(0 < x < 1) is given as 



x)cl 



(3) 



We assume that the initial core is a Bonnor-Ebert sphere, 
whose mass Mq is set to 3 Mq or 30 Mq. A Bonnor-Ebert 
sphere is a self-gravitating isothermal sphere in hydrostatic 
equili brium with a surrounding pressurized medium (IBonnoi 
19561) . Gravitational stability of a Bonnor-Ebert sphere is 
characterized by a non-dimensional physical truncation radius 
^max = rout -sJA-kGpc/c, whcrc ro,„ is the truncation radius, G is 
the gravitational constant, and p^ is the central density. If ^,„„i 
is smaller than the critical value ^crii - 6.45, the Bonnor-Ebert 
sphere is stable. All Bonnor-Ebert spheres which we use in the 
calculations are gravitationally stable {^max - 3,4,5). 

It is reasonable to assume that the initial density distribu- 
tion of the core has spherical symmetry. Density distribution 
of the core is assumed to remain unchanged while the core is 
immersed within an expanding H ii region. The sound crossing- 
time of the Bonnor-Ebert sphere in model C is 7.02 x 10^ yr, 
and the expansion velocity of the H ii region is approximated 
well with Rj = CiiRj/R,,)-^'"^ at Ri > R,„ where Ri and R,, 
are the radii of the H ii region and the Stromgren sphere. If we 
assume /?, = 10 pc and an ambient density of 100 cm"^, the 
timescale on which the Hii region engulfs the core, IR^/Ri, is 
4.92 X 10"* yr. This timescale is smaller than the sound-crossing 
time, and the adoption of spherical symmetry is acceptable for 
simplicity. 

We have computed seven models of collapse of cores ex- 
posed to UV radiation from an O star. We summarize the pa- 
rameters of our simulations in Table[T] For the three models (A, 
B, and C) we assume that mass of the Bonnor-Ebert sphere is 
3 Mq and the ionizing star is an 08V star These three models 
differ in the non-dimensional physical truncation radii of the 
Bonnor-Ebert spheres ^,„a.x- Bonnor-Ebert spheres with larger 
^„M v have higher central densities and smaller radii. For other 
three models (D, E, and F) we assume that the ionizing star 
is an 04V star. In model G, we computed the case of a more 
massive 30 Mq core, exposed to UV radiation from an 08V 
star. 



K. Motoyama et al.: A Radiation Driven Implosion Model for the Enhanced Luminosity of Protostars near HII Regions 



3 



We perform simulations in one spatial dimension using the 
Godunov method. Our hydrodynamical code has second-order 
accuracy in both space and time. The hydrodynamic equations 
are 



Table 1. Parameters adopted in the simulations. 



dp 1 dr pv 
dt r^ dr 

dpv 1 dr^pv^ 



0, 



dr 



dt 

, 2 

dr 



dP 
' dr 



GM(r)p 



(4) 



(5) 



(6) 



where r is the radius, p is the density, v is the radial velocity, 
P is the pressure, and M is the mass within radius r. For near- 
isothermality, the ratio of specific heats, y, is set to be 1.001. 
The gas is assumed to be ideal gas. 

We solved the radiative transfer equation for the diffuse 



UV ra diation using a two-stream approximation. iCanto et al 



(119981) first treated the dif fuse UV radiation usi ng an iterative 
two-stream approximation. lPavlakis et al.l (1200 lb simplified the 
method and showed a non-iterative version serves as a good 
approximation for centrally-peaked density distributions. We 
adopted the latter for our Bonnor-Ebert spheres. The equations 
of evolution for the ionization fraction and propagation of dif- 
fuse UV radiation are given by 



Dx 
ISt 



-ax n + 2(1 - x)crF(r), 



dF 

— = -2(1 - x)no-F(r), 
dr 



(7) 



(8) 



2.7 X 10 '■'cm-' s ' is the recombination coeffi- 

3.0 X 10-"*cm2 



where n is the number density, F is the number flux of ionizing 
photons, a 

cient into the excited state (n > 2), and cr 
is the ionization cross-section. These equations were included 
in the hydrodynamic solver as follows. First, we solved equa- 
tions (|4|i-(|6]) by Godunov method, in which electron density xn 
is treated as a passively advected va riable. After th at, we in- 
tegrated equations (|7])-(^ implicitly ( Williamslll999!) . Pressure 
and internal energy were updated using the new ionization frac- 
tion. 

We ve rified our procedur e abov e by a test problem de- 
scribed in iLefloch & LazareffI (119941) . We simulated the time 
evolution of uniform neutral gas illuminated by an ionizing flux 
increasing linearly with time in one dimensional Cartesian co- 
ordinate. An ionization front propagates with constant veloc- 
ity, and neutral gas accumulated a head of the ionization fron t. 
We adopted same parameters as iLefloch & Lazarefil (119941) . 
Deviations between numerical results and the analytic solutions 
were less than 1 % and 5 % for the density and the velocity, re- 
spectively. Moreover, we checked the accuracy of our code in 
spheric al coordina tes by comparing with the self-similar solu- 
tion oflSh3 (Il977h . 

The sink cell method (e.g.|Boss_& Blacklll982h was used to 
avoid the long computational time encountered near the cen- 
ter. As the gravitational collapse proceeds, the infall velocity 
increases at central region and the time step that satisfies the 
Courant-Friedrichs-Levy (CFL) condition, which is necessary 



Model 


a 


P.' 




Mo' 


spectral type of 


d'' 






[gem ^] 




[Mq] 


ionizing star 


[pc] 


A 


3 


1.05 X 10- 


-20 


3 


08V 


3-10 


B 


4 


3.01 X 10" 


-20 


3 


08V 


3-10 


C 


5 


6.09 X 10- 


-20 


3 


08V 


3-10 


D 


3 


1.05 X 10 


-20 


3 


04V 


3-10 


E 


4 


3.01 X 10- 


-20 


3 


04V 


3-10 


F 


5 


6.09 X 10- 


-20 


3 


04V 


3-10 


G 


5 


6.09 X 10- 


-22 


30 


08V 


3-10 



" Non-dimensional truncation radius of initial core (Bonnor-Ebert 
sphere). 

* Central density of initial core. 
' Total mass of initial core. 
Distance from ionizing star. 

for numerical stability, reduces. If the central density exceeds 
the reference density psM = lOVc, the central 10 grids within 
10 AU will be treated as sink cells. The mass that enters the 
sink cell is treated as a point mass located at r = 0. Boundary 
condition at the outer boundary of the sink cell was given by 
extrapolation from the neighboring cells. Since the infall ve- 
locity near the sink cell is supersonic, this method does affect 
the flow of outer gas. 

We used non-uniform grids in order to obtain high resolu- 
tion at the central region. Sizes of the ith grid from center were 
given as Ar, = (1 -H e)Ar,_i, where e is 0.02, and size of the 
most inner grid is 1 AU. The outer boundary is located at 3rout- 
For example, we allocate 1473 grids in model C. For the core 
itself, 1200 grids were allocated. Reflecting and free boundary 
conditions were imposed at the inner and outer boundaries, re- 
spectively. 

3. Results 

3.1. Time Evolution 

To illustrate common features of the time evolution, we de- 
scribe results in model C with d = 9 pc listed in Table [T] Fig. 
[Ushows the time evolution of density, velocity, ionization frac- 
tion, and pressure. Initially, the core was in hydrostatic equi- 
librium with the rarefied warm neutral gas. At f = 0, UV radi- 
ation, whose UV photon flux is Fj, turned on and illuminated 
from the outer boundary. The outer layer of the core was ion- 
ized by UV radiation. Temperature and pressure at this outer 
layer increased drastically. The hot ionized gas then expanded 
outwards, accelerated to ~ 40 kms"', and drove a Shockwave 
into the core. At f = 1.02 x 10^ yr, the Shockwave reached cen- 
ter of the core. The dashed line in Figdfa) shows that the core 
is compressed into a small region of a radius ~ 3700 AU at this 
time. The density becomes ~ 10^ times higher than the initial 
density and the core becomes gravitationally unstable. 

The core collapses subsequently. We see from FiglT](a) and 
(b) that density and velocity of the compressed core increase 
rapidly after t = 1.02 x lO'' yr due to self-gravity. After the 
central density exceeds the reference value psink, the sink-cell 
treatment turns on, at which point a protostar forms at the cen- 
ter of the core. At this phase, inner region of the core accretes 
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Fig. 1. Time evolutions of (a) density profile, (b) velocity profile, (c) ionization fraction and (d) pressure profile in the model with 
^max - 5 and d - 9 pc. The dotted lines denote the profiles at f = 0. The profiles at f = 0.66 x lO^yr, t - 1.02 x lO^yr, and 
f = 1.09 X lO^yr are labeled by dash-dotted, dashed, and solid lines, respectively. 



onto the central protostar, and the ionized outer layer expands 
outward. 



[3] shows Q as a function of radius. ISpitzeii (|1978h derived an 
analytic estimation of Q as 



Velocities of the shock front and the ionization front remain 
nearly constant, on the other hand. Fig.|2]shows locations of the 
ionization and the shock fronts as functions of time. Velocity of 
the shock front is a constant value of ~ 1.37 kms"'. Analytic 
estimation described in section 14.11 shows a similar velocity 
of 1.18 kms For convenience, we define ionization front 
at the location where the ionization fraction is 0.9. How ion- 
ization front is defined in fact hardly affects the results, since, 
as shown in Fig. \llc), the transition layer from neutral to ion- 
ized gas is very thin. Average velocity of the ionization front is 
1.23 kms"', and analytic estimation shows a similar velocity 
of 1.14 kms"'. 



1 + ji + 



(9) 



where rjf is the radius of the ionization front. Our numerical 
results show a slightly large r value than t hat of the analytic 
estimation by assumptions. Spitz eij ( 1978 ) assumed that the 
photo-evaporated gas expands at the constant sound speed c,. 
As shown in FigHtb), however, the photo-evaporated gas accel- 
erated to ~ 3c, (~ 40k ms"') in our simulations due to pressure 
gradient. This high velocity reduces the density of the photo- 
evaporated gas. As a result, the recombination rate of hydrogen 
atoms is smaller than that of the analytic estimate. 



Some fraction of the incident UV photons is absorbed in 
the layer of the photo-evaporated gas. In this layer, ionized hy- 
drogen atoms recombine with electrons at the rate of ax^n per 
unit volume per unit time. Some UV photons are used for ion- 
ization of these recombined hydrogen atoms. We define Q, the 
ratio of UV photon flux at the ionization front to the incident 
UV photon flux, as the photo-evaporation rate of the core. Fig. 



3.2. Accretion Rates 

As a reference, we simulate the spontaneous collapse of the 
core without UV radiation. Since the thermal pressure balances 
the gravitational force in the Bonnor-Ebert sphere, we slightly 
enhance the density to p = 1.05 pbe, where pbe is the density 
of a marginally stable Bonnor-Ebert sphere {^max = 6.45), to 
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t [10= yr] 

Fig. 2. Locations of the ionization front and shock front. The 
thick (black) and thin (red) solid lines label the position of 
shock front and ionization front in model C with d - 9 pc, 
respectively. The thick (black) and thin (red) dashed lines label 
the analytic estimates described in section 1431 respectively. 




0.6- 

0.4 '- '- 
0.2 - - 

0-0 F • 

10000 20000 30000 

r [AU] 

Fig. 3. The ratio of UV photon flux at ionization front to in- 
cident UV photon flux as a function of radius in the model C 
with d = 9pc. The dashed line indicates the analytic estimation 
which is given by Eq. (|9|l. 



initiate collapse. This core whose mass is 3 Mq collapses by 
self- gravity, and the averaged value of accretion rate is 7.4 x 

10-^ Mq yr"^ 

The RDI indeed increases accretion rates. Fig.|4]shows the 
accretion rate as a function of time in model C with li = 9 pc. 
The accretion rates are measured at the radius of 10 AU. The 
horizontal line indicates the average accretion obtained from 
the reference model. At f = L02 x 10^ yr, the central density 
reaches the threshold density psM, and the sink cell method 
turns on. The accretion rate increases rapidly afterwards, and 
lasts for ~ lO'^yr. Fig.|5]shows the averaged values of accretion 
rates as a function of distance from the ionizing star The accre- 
tion rates increase with decreasing distance from the ionizing 
star, and with the increasing UV photon flux. If the cores are 



10"^ 



T3 




0.9 1.0 1.1 1.2 

t [10= yr] 

Fig. 4. Accretion rate as a function of time in the model C with 
d - 9 pc. The horizontal line indicates the averaged value of 
accretion rate in model without UV radiation. 



located within 10 pc from an ionizing star, the accretion rates 
can go 1-2 orders higher than models without UV radiation. 

The accretion rates not only depend on the UV photon 
fluxes but also on ^„„^ and Mq. Cores with larger values of ^„a.i 
have smaller radii and higher densities. Shock velocities and 
densities in post-shock regions are larger in the models with 
small ^max- Gas of higher density has a shorter free-fall time, 
and reaches an accretion rate higher in the models with smaller 
^mov We will compare this with analytic estimates for better 
understanding. In the model G, not only small free-fall time 
but also large mass of compressed core is responsible for high 
accretion rate. 

3.3. Final Masses 

Photo-evaporation can reduce the actual mass that accretes 
onto the protostar. In the outer layer, hot ionized gas expands 
outward at a velocity up to ~ 40 kms ^ Photo-evaporation 
of gas takes place at s/lGM/r, except for the case of very 
small cores. The escape velocity of the Bonnor-Ebert sphere 
in model C is ~ 0.42 km s \ for example, is much smaller than 
that of the photo-evaporated flow. Photo-evaporation prevents 
large amount of gas from accreting onto the protostar, and final 
mass of the protostar can be small due to mass loss by photo- 
evaporation. 

Final masses of protostars decrease with decreasing dis- 
tance from the ionizing star Fig. |6] shows the final masses 
of protostars as functions of distance from ionizing star The 
mass of a protostar is defined as mass which has entered sink 
cells, and calculations continue until accretion rates become 
sufficiently small (10^^ MQyr"'). Mass loss due to photo- 
evaporation at the surface of the core is proportional to the inci- 
dent UV photon flux and final masses of protostars are smaller 
near the ionizing star. For an 08V or an 04V star, mass loss 
due to photo-evaporation varies from several percent to several 
tenth of percent of the initial core mass. The more massive the 
ionizing star gets, the smaller the final mass. 
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d [pc] 




2 



4 6 8 
d [pc] 



10 12 



Fig. 5. The averaged accretion rates as functions of distance 
from the ionizing star Open diamonds, open squares, open cir- 
cles label results from the model A, B, and C, respectively. 
Filled diamonds, filled squares, filled circles label models D, 
E, and F, respectively. The crosses are for the model G. The 
horizontal line indicates the averaged accretion rate in the ref- 
erence model. Dotted line, dash-dotted line and dashed line in- 
dicate analytic estimation for models A, B, and C, respectively. 



Fig. 6. Ratio of final mass of protostar to initial core mass as 
a function of distance from ionizing star The open diamonds, 
the open squares, the open circles indicate results in the mod- 
els A, B, and C, respectively. The filled diamonds, the filled 
squares, the filled circles indicate results in the models D, E, 
and F, respectively. The crosses indicate results in the model G. 
The dotted line, the dash-dotted line and the dashed line indi- 
cate analytic estimation for model A, B,and C, respectively. 



The final masses not only depend on the UV photon fluxes 
but also on the sizes of cores. Fig.|6]shows that mass loss due to 
photo-evaporation is significant in the model G, in which initial 
radius of the core is 10 times larger than that in the model C. 
Larger core has larger photo-evaporation rate and longer time 
to be exposed to UV radiation because of larger surface area 
and longer crossing time of shock wave. This is the reason why 
the core lose large amount of mass in the model G. For the same 
reason, final masses are small in the models with small ^max- 

4. Discussion 

4.1. Analytic Estimates 

In this section, we derive analytic estimates of final mass and 
accretion rates of protostars under the influence of RDI. For 
simplicity, we assume that a core has a uniform density no and 
a radius Ro initially. The core is exposed to UV radiation whose 
number flux of UV photons is F,. Spherical symmetry is as- 
sumed. 

The final mass of protostar Mj- is expressed as 



Mf = M( 



■ 



4nrif(t) mdt. 



(10) 



where rif(t) is the radius of the ionization front at time t, m 
is the mass loss flux at radius rifit), and is the crossing 
time of the Shockwave. The second term on the right-hand side 
of Eq. ( fTOl ) expresses mass loss due to photo-evaporation. We 
only consider mass loss during the compression phase when 
the core is compressed into a small region by the shock. After 



the compression, the core collapses within the free-fall time of 
the compressed core that is much shorter than f,. We need m, 
ts, and rjf{t) to estimate the final mass of protostars using Eq. 

Coll. 

The mass loss flux ih due to photo-evaporation is expressed 

as 



m = ninQFi. 



(11) 



where Q is the ratio of UV photon flux which arrives at surface 
of core to incident UV photon flux F,. The UV photons ionize 
the outer layer of the core, but some fraction gets absorbed in 
the photo-evaporated flow. As shown in Fig. [3l Q approaches 
to unity as the ionization front enters the inner region of the 
core. However, we assume that Q remains constant through the 
compression phase and calculate the value of Q at the radius R{) 
using Eq. Q. 

The shock velocity m , is derived following 'Bertoldi (Il989h 
to estimate the shock crossing time fj. Assuming ram pressure 
of the flow equals thermal pressure of the post-shock region, 
we obtain 



(12) 



where no and ni are number densities of the undisturbed region 
and the post-shock region, respectively. The equation of mass 
continuity at the ionization front gives 



QFi = niu',p. 



(13) 



where u'j^ is the velocity of ionization front in the rest frame of 
the post-shock region. The ionization front is assumed to be D- 
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critical (cf. lSpitzeij|l978b . and velocity of the iotiizatioti frotit 4.2. Luminosity 



IS given as 

^2 



U 



IF 



2c,- 



Substituting Eq. ( fT4b into Eq. ( fT3] ) gives 
2QciFi - clni. 

Combining Eq. ( fTSI ) and Eq. ( fTST i, we obtain 



2ciQFi 
no 



Using Eq. ( fT6b , the shock crossing time is expressed as 



V 



no 



(14) 



(15) 



(16) 



(17) 



We derive the velocity of ionization front ujf in the rest 
frame of the undisturbed region to obtain radius of the ion- 
ization front rjf(t). Taking into account ni - «o and 
riQUs - ni(Us - vi), we obtain 



Vl - Us 



(18) 



where vi is velocity of the post-shock region in the rest frame 
of the undisturbed region. Velocity of the ionization front is 



2 2 

c c 

U,F = M.V + — . 

Us 2Ci 



(19) 



Using u/f, radius of the ionization front is expressed as 
riF - Ro - uift. (20) 
The final mass of protostar using Eq. ( fTOl l gives 

u^ 

Mf = Mo - AnmnQFi (Rlts - Rou,Ft] + (21) 

This equation gives us the analytic estimate of final mass of a 
protostar influenced by RDI. 

The accretion rate, over a free-fall time, is roughly esti- 
mated by 



Mf 
M = — , 



(22) 



where /// - -^5^^ is the free-fall time of the compressed 
core. The core is compressed into a radius of rjf(ts) when the 
shock reaches its center. The density of the compressed core is 
then given as 



Mf 



(23) 



Luminosity of a protostar is mainly produced by accretion, be- 
cause nuclear burning does not influence much of the luminos- 
ity during the protostellar collapse phase. The accretion lumi- 
nosity L is expressed as 



GM,M 



(24) 



where M, is the mass, M is the accretion rate, and Rt is the 
radius of the protostar Enhancement of accretion rates due to 
RDI has a large effect on the luminosity. Although the accretion 
rate depends not only on the distance from the ionizing star 
but also on the ^^ax of the Bonnor-Ebert sphere, it is relatively 
not as sensitive to ^^ax compared with the distance. Our results 
show that for cores located within 10 pc of an ionizing star, the 
accretion rates becomes 1-2 orders of magnitude higher than 
those without triggering from the UV radiation. 

Enhanced luminos ity around bright-ri mmed clouds is 
seen in observations. ISugitani et al.l (11989 ) observed three 
bright-rimmed globules associated with cold IRAS sources in 
'^CO(J=1-0) emission. They compared the physical parame- 
ters with four isolated dark globules and found that the ra- 
tios of the luminosity to globule mass, 3-13 Lq/Mq, are 
1-2 orders higher than the values of isolated dark globules, 
0.03 - 0.3 Lq/Mq, as shown by our calculations. These ob- 
jects are good candidates influenced by the RDI. 

Protostars near H n regions may spend most of their lifetime 
in the state of high luminosity. If the duration of high accretion 
rate is much shorter than the lifetime of the protostar, it is rela- 
tively diflicult to catch the protostar in its high-luminosity state 
by observations. In the simulations, the protostellar phase be- 
gins from the time when the sink-cell method is turned on and 
lasts until 99 % of the final mass has accreted within the sink 
cells. The lifetime of the protostar in model C is 9.1 x 10^ yr 
with c/ = 9 pc. As shown in Fig. 21 the accretion rate is en- 
hanced for a duration of 8.4 x 10^ yr, comparable to the lifetime 
of the protostar. It is likely that protostars near H 11 regions be 
observed in the state of higher luminosity. 

4.3. Comparison with observations 

Evidences of triggered star formation have been suggested on 
edges of many H 11 regions. Some of the observations in fact fa- 
vor the scenar i o of s equential star formation around a massive 
star. lLee et al J (120051 ) selected candidates of pre-main sequence 



stars based on 2MASS colors in the Orion region, and revealed 
spatial distribution of classical T Tauri stars using their spec- 
troscopic observations. The pre-main sequence stars form be- 
tween the ionizing star and the bright-rimmed clouds, but not 
deeply inside the cloud. The youngest stars are found near the 
cloud surfaces, forming an apparent age gradient. As the Hn 
region expands, it might be able to trigger star formation se- 
quentially from near to far Several of such examples are sum- 
marized in Table 2. 

RDI may contribute to the triggering mechanism around a 
massive star. iLee et al.l (l2005b found that only bright-rimmed 



and Eq. {22\ gives us the analytic estimate of the accretion rate, clouds with strong IRAS 100 fim and Ha emission show signs 
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Table 2. Star formation regions of potential triggered origins 



Cloud 


Spectral type of 


Projected distance from 


Estimated shock 


Evidence of triggered 




main ionizing star 


ionizing star 


velocity 


star formation 


molecular pillars in Ml 6' 


05V 


2 pc 


1.3 kms-' 


age gradient of YSOs 


IC1396N- 


06.5f 


11 pc 


0.6 kms' 


age gradient of YSOs 


SFO 1 1 , SFO 1 1 NE, SFO llE^ 


06V 


11 pc 


1.4 kms"' 


on-going star formation 










near ionization front 


Orion bright-rimmed clouds'' 


08IIF or Oe:* 






on-going star formation 










near ionization front 



" Ionizing star of B30 and B35. 

Ionizing star of IC2 1 1 8, LDN 1 6 1 6 and LDN 1 6 34. 

References: (l) lFukuda et"d] ( l2002h . lwhite et al.l ( ll999l) . (2) lGetman etal] ilOOlt) . (3') lThomDson et alj ( l2004l) . (4) lLee et al.l ( l2005h 



of ongoing star formation. This suggests that star formation 
may have been triggered by strong compression from the ion- 
ization shock, because these emissions are tracers of dense gas 
and ionization. Small-scale age gradients of YSO s near the ion 



ization fronts are seen in some molecular clouds. Fukuda et al 



(12002') observed heads of molecular pillars Hi and 112 in the 
Eagle Nebula (M16). They found that protostars and starless 
cores are ordered in age sequences. The more evolved objects 
are closer to the cloud heads (and the ionizing star) and the 
younger objects are farther from the exciting st ar. Similar age 
gradi ents are also seen in bright-rimmed clouds (jSugitaiii et al 



19951) and cometary globule 1C1396N (Getm anet all 120071) . 
These age gradients near ionization front seem to suggest that 
star formation may be sequentially triggered by ionization 
shocks. 

Stars formed facing th e H ii regions appear to be more lu 



Yamaguchi etal] ^999) carried out '^CO(J =1-0) 



minous. 

observations toward 23 southern H n regions and identified 95 
molecular clouds associated with Hit regions. They found that 
IRAS point sources tend to be more luminous on the side fac- 
ing an H II region. Their result is consistent with enhanced lu- 
minosity as a result of the RDl. 

Shock speeds derived from observations are similar. 



White et all (11999!) and [Thompson etal.1 (120043 derived shock 
velocities from pressure in pre- and post-shock regions. 



Getman et al.l (l2007h derived shock velocity from the age gradi- 



ents of YSOs. Their estimated values of a few km s are con- 



sistent with our results. For example, [Thompson et al.l (120041) 
estimated that shock velocity in bright-rimmed cloud SFOll 
is ~ 1.4 kms SFOll is illuminated by an 06V star whose 
projected distance is 1 1 pc. These conditions are close to those 
in model C with d - 9, whose ionizing star is an 08V star 
As described in section 13.11 the averaged shock velocity is 
~ 1.37 km s ' in this model. 



Lee & Chen (l2006l) found that classical T Tauri stars ex- 



ist near the surfaces of bright-rimmed clouds and intermedi- 
ate mass stars exist relatively inside clouds. One possible dif- 
ference of stellar mass may be density or mass in the parent 
cores. Photo-evaporation may reduce the mass that would oth- 
erwise accrete onto the stars. In our simulations a few percent 
to several tenth percent of initial mass may be lost by photo- 



evaporation. Mass of young stars may be smaller near the sur- 
face of bright-rimmed clouds compared with inner regions. 

5. Conclusions 

We simulate the evolution of cores exposed to strong UV ra- 
diation around the Hit regions. We investigate possible effects 
of the radiation-driven implosion on the accretion luminosity 
of forming protostars. We also developed analytic estimates of 
accretion rates and final masses of protostars, which are com- 
pared with the numerical results. 
The main findings are: 

1 . Radiation-driven implosion can compress the cores and de- 
crease the free-fall time. 

2. The accretion rates are positively influenced by on the inci- 
dent UV photon flux. 

3. Photo-evaporation of the parent cores may decrease the fi- 
nal masses of protostars. However, mass loss due to photo- 
evaporation does not affect significantly the accretion lumi- 
nosity. 

4. Final mass of a protostar depends mainly on the size of 
the parent core and incident UV photon flux. Final mass of 
a protostar decreases with the increase of UV photon flux. 
Dense compact cores, on the other hand, are hardly affected 
by mass loss due to RDl. 

5. The RDl increases luminosity of protostars by orders of 
magnitude than spontaneous star formations without RDl. 
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